Soliton Fermi Sea in Models of Ising-coupled Kondo impurities 
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We study a model of Ising-coupled Kondo impurities that can be applied to quantum dots with 
capacitance coupling, coupled qubits with an incoherent environment, etc. We show that this model 
becomes equivalent to a Anderson impurity model with a novel solitonic Fermi sea. We derive exact 
results at the Toulouse point, and use the flow equation method to extend this analysis away from 
the Toulouse point for not too large Ising couplings. 



Introduction. Kondo physics plays a fundamental 
role for the low-temperature behavior of a large vari- 
ety of physical systems like magnetic impurities in met- 
als, heavy-fermion systems, glasses, quantum dots, etc. 
Its key feature is the quenching of the impurity entropy 
through nonperturbative screening by many-particle ex- 
citations in the associated quantum bath. For magnetic 
impurities in metals this amounts to the formation of the 
Kondo singlet between the localized spin and electron- 
hole excitations in the Fermi sea^. 

Most of the aspects of single-impurity Kondo physics 
are now well understood after theoretical tools have been 
developed that can deal with its intrinsic strong-coupling 
nature^. However, in many physical systems the inter- 
action of different impurities, i.e., multi- impurity Kondo 
physics, is important. For example in heavy-fermion sys- 
tems the RKKY interaction between different impurity 
spins is responsible for the competition between local 
Kondo physics and long-range magnetic order that de- 
termines their phase diagram^, or in glasses the phonon- 
mediated coupling between different tunneling centers is 
known to be importantS^. More recently, related ques- 
tions about coupled two-level systems have gained much 
interest in quantum computation, where decoherence due 
to unwanted couplings among qubits and between qubits 
and environment should be avoided; on the other hand 
the intentional coupling of qubits is the key step to per- 
forming quantum logic operations. 

In the present paper we investigate the case of two 
Kondo impurities coupled via an Sz-Sz interaction, i.e., 
an Ising-type coupling. Equivalently, one can think of 
two two-level systems with transversal coupling, with the 
experimental realization of coupled flux qubits^. Based 
on the flow equation method^ we can derive analytical 
results for this problem by mapping it to a conventional 
Anderson impurity model. This mapping allows us to 
derive exact results at the Toulouse point, and controlled 
approximations in its vicinity. Also we will show how our 
results can be extended to an arbitrary number of Ising- 
coupled impurities with an inflnite-range interaction. 

Coupled impurities or two-level systems have been in- 
vestigated in a number of papers2ii2iiiii2ii^, however, 
most attention has been focussed on the case of SU(2)- 
symmetric coupling between the impurity spins. We will 
show that new phenomena arise in the Ising-coupled case 
due to the two-fold degeneracy of the ground state of the 



isolated impurity systemic. In the context of Kondo im- 
purities, this Ising-like coupling can be thought of as an 
effective impurity interaction for heavy-fermion systems 
with an easy axis. Also, Ising coupling appears naturally 
in quantum dots that are coupled via their mutual capac- 
itance: here the two-level systems are pseudospins repre- 
senting the number of electrons on the dots, and therefore 
SU(2) symmetry is broken from the outsetiS*iSJi. 

Model. Let us first formulate the model in the lan- 
guage of Kondo impurities. For simplicity, we will con- 
sider a situation where the two impurities are placed 
in different fermionic baths, yielding the Hamiltonian 

= Hf + H 
ual Kondo Hamiltonians 



H^2^^- where HY" are the Individ- 
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with the fermion creation/annihilation operators 



(1) 



''koti ' 



^kai with wave vector k 
and belonging to Fermi sea 



2'Kn/L, n G Z, spin index a 



denote the 

localized electron state at the impurity sites and Si, S2 
are the two spin-1/2 impurities. We will assume a linear 
dispersion ek — v k throughout this paper and set ep = 0. 
The Ising coupling, which breaks SU(2) symmetry, takes 
the form 



H 



Ising 
12 



K SI 



(2) 



Note that i?j2 ^ doubly degenerate ground state for 

both signs of K . The simplifying assumption of two dis- 
connected baths is naturally satisfied for coupled quan- 
tum dots; for magnetic impurities in metals it misses the 
RKKY interaction, which we instead explicitly include 
by H12. 

The case of SU(2)-symmetric inter-impurity coupling, 
K Si ■ S2, has been the subject of numerous investiga- 
tions in the pasliSiiSiiiiiSiH. Generically, two different 
regimes are possible as function of the inter-impurity ex- 
change K: for large antiferromagnetic K the impurities 
combine to a singlet, and the interaction with the con- 
duction band is weak, whereas for ferromagnetic K the 
impurity spins add up and are Kondo-screcned by con- 
duction electrons in the low-temperature limit. In both 
parameter regimes the ground-state spin will be zero. 
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and there is no quantum phase transition as K is var- 
ied in the generic situation without particle-hole sym- 
metry (whereas one finds an unstable non-Fermi liquid 
fixed point in the particle-hole symmetric case)i24ii*i^. 
However, the finite-T crossover behavior is significantly 
different in the two regimes: for large antiferromagnetic 
K, the singlet is formed at an energy scale of order K. 
In contrast, for ferromagnetic where both spins have 
to be quenched, a two-stage Kondo effectSii^ occurs if 
the two screening channels have different Kondo coupling 
leading to different Kondo temperatures: upon lowering 
the temperature the two spin- 1/2 degrees of freedom are 
"frozen" at two distinct temperature scales. 

An alternative starting point for our model of Ising- 
coupled impurities can be formulated in terms of two 
two- level systems (spin-boson models), H^^ = Hf^ + 

^SB 1 u^trans 



Tjtr 



with 



H: 



SB 



A,: 



k>0 



(3) 



fe>0 



and transversal coupling between them, Hl™'^ = 
^ erf erf. Here b\^, 6^,^ are the bosonic cre- 
ation/annihilation operators for heat bath #i. is the 
bare tunneling matrix element for the two-level system 
described by cr^, and hi is the asymmetry of its two wells. 
The impurity properties are completely parametrized by 

the spectral function Ji{uj) A^^ 5{uj — ujki), which 

we assume to be of Ohmic form, Ji(w) = 2aiUJ e'^^^"" . 
One realization of this model is the interaction of tun- 
neling centers in glasses through higher-order phonon 
exchange^. In the context of quantum computation 
this model arises in studies of decoherence of coupled 
solid state qubits: the transversal coupling K is gener- 
ated through a superconducting flux transporter, and 
the heat baths describe the environment leading to 
decoherence^'^. 

With a polaron transformation H^^ can be mapped 
onto the following Hamiltonian with the same impurity 
properties, H^°^ = H^°^ 



jjhos 



^trans^ Here 



- I hot + ft [VdK; 0) ar + h.c.) + ^ vkblfi^ 



k>0 



ki 



(4) 



with coupling constants ft = Ai/-v/8 and ujki — vk with- 
out loss of generahty. T^(A; a;) denote vertex operators 



V, 



(A; .) exp [aJ: e-^'/\e^^'^%^ - e'^^^^.J 



fc>0 



with scaling dimension Ai — ^/2al and a = v/ojc- 

Eq. is also the bosonized form describing the low- 
energy properties of the Kondo Hamiltonian ^ with 
^L' ^fei describing the spin-density excitations in Fermi 
sea #i, a = ppv, Xi — \/2 — {pfJ)/VSt^, } 



(ji = J/2. In the sequel we will use Q as the unified 
starting point of our analysis and use the terminology of 
Kondo impurity problems. Notice that H^^^"^^ remained 
unaffected by the above unitary transformation. 

Toulouse point. In a first step to study the phase di- 
agram of Ising-coupled Kondo impurities we investigate 
the special Toulouse point^* Ai = 1 [ai = 1/2), where 
we will show that the problem can be solved exactly. 
Since the vertex operators ^^(±1; x) obey fermionic anti- 
commutation relations, the Hamiltonian of an individual 
Kondo impurity, Hf°^, can be refermionized, leading to 
a resonant level model of spinless fermions for each of the 
two impuritiesiS. The Ising coupling transforms into a 
density-density interaction, and the full model H^°^ thus 
reads 

= ^ « fc ele,^ + Vk^ {eld, + dfe.J 



ki 



^edindi + nd2 - 1) + {udi - nd2) 



(5) 



with the representation a^' — d|, a., = , af = 2{dld^ 
1/2), rid, = Ffc 
hoc ^ hi - h2, and 



giy/2TTa/L, = {hi + /i2)/2, 



f dcf 



1 



dx e^"" 



y^{\■.x) 



(6) 



0, and 



Eq. (0) represents one of the main results of our paper: 
At the Toulouse point the problem of two Ising-coupled 
Kondo impurities maps onto an Anderson impurity model 
(AIM) with Coulomb interaction K. The Fermi sea of 
this effective AIM consists of the soliton excitations © 
constructed from spin-density excitations around the im- 
purities. The index i — 1,2, originally labelling the two 
impurities and their baths, now plays the role of a pseu- 
dospin in this effective AIM. The impurity properties 
are completely described by the hybridization function 

Ai(e) =^ J2k ^ki <^(e— £fe); which is constant in the above 
model Ai(e) = pF9i- With this mapping the original 
problem (at the Toulouse point) is essentially solved since 
(O is well-understood using various methods like, e.g., 
Wilson's numerical renormalization group^, and for the 
special case of both impurity sites being equal {Vki = V 
and hioc = 0) is even integrable with exact Bethe ansatz 
results^. 

This integrable case allows us to study in detail the 
competition between local Kondo screening and the Ising 
interaction. From the Bethe ansatz solution one knows 
that the ground state of the Anderson impurity model 
is always unique with no phase transition as a function 
of K. For simplicity we will first analyze the problem 
with no local magnetic fields hi — h2 ~ (cd = 0). For 
large positive K ^ A(0) (large antiferromagnetic Ising 
interaction) one finds (n^i -|- 71^2) 1 {{Sf + 51) 0): 
the impurity site of the effective Anderson model is singly 
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Temperature T 



FIG. 1: Schematic temperature dependence of the total im- 
purity entropy of two Ising-coupled Kondo impurities for the 
case of large inter-impurity coupling, K ^ T^. Below K, the 
impurity spins combine into the two possible (anti)parallel 
low-energy states for (anti)ferromagnetic K. The remaining 
pseudospin degree of freedom is screened below T]^". 



occupied. There are two degenerate impurity configura- 
tions, namely ti and It- However, these two configura- 
tions get mixed by many-particle excitations in the effec- 
tive Fermi sea. Thus, the pseudospin is screened below 
a new collective Kondo scale T^^^ ~ exp(— |i4r|/8A(0)), 
and the ground state is a singlet. On reducing the tem- 
perature the impurity entropy is quenched in two stages 
(Fig. Pl: for r > A' one finds S'imp = In 4, which is 
reduced to the two singlet configurations (S'imp = In 2) 
for temperatures of order K, and finally quenched com- 
pletely below the collective Kondo temperature T^^^. 
Note that this type of two-stage screening is completely 
different from the one mentioned in the introduction^, 
which occurs if the two conventional Kondo screening 
channels have different strengths. 

For large negative K (large ferromagnetic Ising interac- 
tion) the situation is related to the antiferromagnetic case 
by particle-hole symmetry: one finds K-S'i+S'l)! 1, i-e-, 
the two impurity states being degenerate in the absence 
of the baths are It and J,|. The screening of the asso- 
ciated pseudospin appears at the same collective Kondo 
scale as above. 

For smaller \K\, namely of order A(0) or below, the 
only relevant temperature scale is set by A(0): the two 
impurities do not couple into an inter-impurity doublet 
since they are screened individually by their respective 
Kondo interaction at temperatures above K. Notice 
that at the Toulouse point A(0) can be identified with 
the Kondo temperature of the individual impurity i, 
below which screening occurs. 

Going away from the integrable case the effective 
Anderson impurity model Q still displays the same 
crossover behavior between the various limiting cases as 
a function of K, provided that K fulfills Eq. H1U|) below. 
In general the above crossover temperatures then depend 
on Cd and the local magnetic field ft-ioc. 

Flow equation solution. Next we generalize the exact 
results at the Toulouse point to general couplings using 



the flow equation diagonalization of the Kondo Hamilto- 
nianSS. In Ref.|22| a continuous sequence of unitary trans- 
formations was constructed that led to a systematic ap- 
proximation scheme which becomes exact at the Toulouse 
point and yields a controlled approximation away from it. 
This means we know a unitary transformation Ui that di- 
agonalizes the Hamiltonian Q up to higher-order terms 
that do not affect the universal properties. It was shown 
that Hf""" = [/. Hf""" U} takes the form 

fc>0 

+ Y,ujk [V{Uk)-k),V{-Uk);k)] , 

k 

where the third term consists of Fourier transformed ver- 
tex operators with a scaling dimension \i{k) that de- 
pends on the energy scale. The unitarily transformed 
operator = cr^Uj is given bjsSS 

a,f =^afeafc4y(A,(fc);fc),y(-A,(fc');fc')] (8) 

k,k' 

with certain coefficients ak- We now make use of this 
single-impurity flow equation diagonalization by applying 
the combined unitary transformation U = U1U2 on the 
coupled system 

H = U.U^HUM (9) 
K 

rrbos , rrbos , x.'^ 

- J^l + J^2 + ^ <^1<^2 ■ 

This procedure is summarized in Fig. [3 In the sequel 
we will only be interested in the strong-coupling regime 
of the single-impurity problems Af < 2 implying a flow 

to the Toulouse point in the infrared limit Ai(fc) — > 1. 
We now approximate JHJ by setting Xi{k) — 1 Vfc, i.e. 
by using fermions instead of the vertex operators with a 
running scaling dimension. This approximation becomes 
exact in the low-energy limit, where many-particle pseu- 
dospin excitations are most important. 

Using this approximation (jSJ can be refermionized 
leading to a Hamiltonian that has the same structure 
like at the Toulouse point, except with a non-constant 
hybridization function Ai(e) that depends on the uni- 
tary transformation Ui. One can show Ai(0) ~ T^, 

and Ai(e) ~ ppgf\pF€\^'~^ for |e| » T^, with a smooth 
crossover in betweenSi. Here is the individual Kondo 
scale of one impurity without Ising coupling to the other 
impurity. 

The above analysis neglects the feedback of one im- 
purity on the other. The accuracy of the approxima- 
tion within the flow equation transformation Ui decreases 
away from the Toulouse point. Therefore the error in 
the transformation of is proportional to the distance 
to the Toulouse point (A — 1), leading to an error in the 
transformation of the interaction between the two Kondo 
impurities that is proportional to jAT (A — 1)|. When this 
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Two Kondo impurities 



^ bosonization, flow eq. transf., refermionization 
Two resonant-level models | 



Ising coupling — > density-density interaction 
Effective Anderson impurity model | 



FIG. 2: Steps of our approximate solution of the two- 
impurity Kondo problem with Ising coupling. Each individ- 
ual Kondo problem is transformed into an effective spinless 
resonant-level model by means of a flow-equation transforma- 
tion. The original Ising coupling maps onto a density-density 
interaction, leading finally to an Anderson impurity model. 



For simplicity we will only look at the Toulouse point 
(also hi — V«), though the above flow equation analy- 
sis of the Kondo regime can be carried over as well. After 
refermionization one finds an A'^-channel Anderson impu- 
rity model describing the SU(iV) pseudospin fluctuations: 



ki ki 

+K ^ (n,, - (n,, - I) . (12) 



energy scale becomes comparable to the energy scale of 
the individual Kondo Hamiltonians our approximations 
are no longer reliable, which means that we can only 
study the phase diagram for Ising couplings with 



(10) 



For not too large Ising couplings K that fulfill this condi- 
tion we can trust the mapping to the Anderson impurity 
model with Coulomb interaction U ^ K, and carry over 
our analysis from the Toulouse point. 

For larger couplings K the fiow equation transforma- 
tion needs to be extended to higher orders. A careful 
analysis shows that a quantum phase transition can oc- 
cur when the condition (|10|l is violated^''. The Anderson 
model becomes supplemented by a density-density inter- 
action between the impurity orbital and the first conduc- 
tion electron site. This interaction maps onto a ferro- 
magnetic spin-spin interaction when charge fluctuations 
are frozen out for large K; this ferromagnetic interaction 
can in turn destroy the screening of the pseudospin and 
lead to a phase with a doubly degenerate ground state. 
For a full account of this analysis and other approaches 
to the model of Ising-coupled Kondo impurities we refer 
to Ref. m 

A'^ impurities. The above mappings can be general- 
ized to the case of N impurities with infinite-range Ising 
interaction 



N 



1=1 



(11) 



Again the crossover between ferromagnetic and anti- 
ferromagnetic coupling K is known to be smooth. For 
the special case of large positive K ~ this situation corre- 
sponds to a fully frustrated Ising model coupled to dissi- 
pative baths - one can perform a generalized Schrieffer- 
Wolff transformation to an SU(A^) Kondo model that 
describes the quenching of the pseudospin fiuctuations 
below a collective Kondo scale T^^^. Away from the 
Toulouse point we find the same picture as long as the 
condition (|10|l is fulfilled. 



Conclusions. Summing up, we have investigated the 
case of Ising-coupled Kondo impurities. At the Toulouse 
point we could map this model exactly to an Ander- 
son impurity model with a Fermi sea that consists 
of fermionic soliton excitations. The Ising coupling 
leads to two degenerate ground states of the isolated 
double-impurity system; this pseudospin degree of free- 
dom is screened below a collective Kondo scale T^°'' 
through soliton excitations. This piece of physics describ- 
ing many-particle excitations of solitons (many-many- 
particle excitations) is absent for an SU(2)-synimetric 
inter-impurity coupling. Using the novel fiow equa- 
tion methoc&2S we could extend this analysis from the 
Toulouse point for not too large Ising couplings K that 
obey Eq. ifTUI) . 
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